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Abstract. Recently, a new type of symmetry for three-qubit quantum states was 
introduced, the so-called Greenberger-Horne-Zeilingcr (GHZ) symmetry. It includes 
the operations which leave the three-qubit standard GHZ state unchanged. This 
symmetry is powerful as it yields families of mixed states that are, on the one hand, 
complex enough from the physics point of view and, on the other hand, simple enough 
mathematically so that their properties can be characterized analytically. We show 
that by using the properties of GHZ-symmetric states it is straightforward to derive 
optimal witnesses for three-qubit entanglement. 
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1. Introduction 

Entanglement is an important resource in quantum information theory. Therefore, the 
possibihty to detect whether or not a given state is entangled is of significant importance. 
For multipartite systems, there are inequivalent types of entanglement. A given task 
generally requires a certain type (or class) of entanglement, and therefore also class- 
specific detection of entanglement is important. 

A common scheme to classify multi-party quantum states is to consider their 
equivalence with respect to stochastic local operations and classical communication 
(SLOCC) [T, The simplest multipartite case with several entanglement classes is 
a system consisting of three qubits for which Diir et al. have identified six different 
SLOCC classes of pure states, and especially two different classes of "true" tripartite 
entanglement, the Greenberger-Horne-Zeilinger (GHZ) class and the W class. The 
representative states of these classes are the standard GHZ state 

|GHZ) = i=(|000) + |lll)) (1) 
and the W state 

\W) = -L(|001) + |010) + |100)) (2) 

which both are invariant under qubit permutations. 

This classification can be extended to mixed states |3j. It reveals that there is a 
linear hierarchy of entanglement classes. In this hierarchy, the GHZ class represents the 
"highest" entanglement type, followed by the W class, the biseparable states B, and 
finally the fully separable states S so that we have GHZ D D i? D s|3 Mixing states 
of a higher class can result in a state of a lower class, but not vice versa. An obvious 
example is that of equally mixing the elements of an orthogonal basis of GHZ states 
(GHZ type), which gives the completely mixed state (fully separable). Note that, by 
definition, each set of states which has at most a given entanglement type is always 
convex. 

The standard method to detect the SLOCC class of a given state are entanglement 
witnesses [21 lU El E] , see also [7] and references therein. Despite the impressive progress 
in entanglement theory during the past decade [H HI [10], to date there is still no 
method which could determine the SLOCC class of an arbitrary three-qubit state with 
certainty. Therefore, exact results for special families of states are important towards 
a better insight into the structure of the space of mixed quantum states. One way of 
obtaining such results is the reduction of the number of free parameters by introducing 
certain symmetry requirements [HI [121 [131 [H]- Recently, a complete characterization 
of the entanglement type in the family of three-qubit GHZ-symmetric states has been 

I In the following we will use the names of the SLOCC classes in an exclusive sense, that is, we will 
refer to the set GHZ\W as GHZ-class states, to the set W \ B as W^-class states, and to i? \ S as the 
biseparable states. 
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achieved pLSj. That family contains in particular the important case of pseudo-pure 
GHZ states, that is, mixtures 

= p |GHZ)(GHZ| + (l-p) h (3) 

of the GHZ state ([T]) and the state Pmixed = moreover < p < 1. The state Pmixed 
is fully unpolarized and is often regarded as a model for white noise. 

However, in experiments one cannot always assume white noise, and for more 
general admixtures the resulting state will not be GHZ symmetric. Therefore, methods 
are required to determine the entanglement class of such states. A possible strategy is 
to select certain subsets of mixed states which lead to a more adequate representation 
of a given type of noise and allows for a more precise determination of an experiment- 
specific optimal witness. This is the approach that we pursue in this article. To this end, 
we first review the properties of GHZ-symmetric three-qubit states. Subsequently, we 
explain the usefulness of witnesses whose optimality is restricted to certain subsets |16j . 
Finally we explicitly demonstrate how GHZ-symmetric states may serve to find better 
witnesses for pseudo-pure states, in particular for the mixtures in equation ([3]). 



2. GHZ-symmetric three-qubit states 

In this section, we summarize the properties of GHZ-symmetric three-qubit states |15j . 
In addition, we determine the set of GHZ-symmetric states with positive partial 
transpose. 

Density matrices with GHZ symmetry are invariant under unitary transformations, 
denoted GHZ symmetry operations, which leave the standard GHZ state equation ([1]) 
unchanged: 

(i) arbitrary qubit permutations, 

(a) simultaneous three-qubit flips (i.e., application of ^ <Jx ^ cTx), 
(Hi) correlated qubit rotations about the z axis of the form 

Here, and o"^ are Pauli operators. Note that the operations (i)-(iii) do not change 
the entanglement of an arbitrary three-qubit state. 

The GHZ-symmetric three-qubit states can be parametrized as 

p{.x, y) = + 7rGHZ+ + - ttghz^ ^ ~ v^) 
where tt^ = is the projector onto state \ip) and we have defined the GHZ-t states 

|GHZ±) = i=(|000)±|lll)) . (6) 

The y coordinate ranges in the interval -1/(473) < y < while, for given y, 

the X coordinate may assume values |x| < (\/3y/2 + 1/8). The states thus form a 
plane triangle (shown in Fig. 1) with corners {x,y) = (— 1/2, -\/3/4) corresponding to 
GHZ_, {x,y) = (+1/2,^3/4) corresponding to GHZ+, and {x,y) = (0, -l/(4v/3)) 
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which represents the state pr = J2i=ooi "^^^ completely mixed state is located at 

the origin {x,y) = (0,0). 




Figure 1. The set of GHZ-symmetric three-qubit states [15]. We have indicated the 
four SLOCC classes of three-qubit entanglement: GHZ (grey), W (yellow), biseparable 
(green), and separable (blue). At the origin, there is the state pmixod — The solid 
magenta line illustrates the pseudo-pure states ([3]). 



The GHZ-symmetric states provide an example to observe the hierarchy of 
entanglement classes mentioned in the introduction. A GHZ-symmetric state is fully 
separable iff it is in the polygon defined by the four corner points (0, — l/(4v^)), 
(1/8,0), (0,4/^3) and (-1/8,0) (the blue area in Fig. 1). Correspondingly, it is at 
most biseparable iff it is in the polygon with corners (0, — l/(4v^)), (1/4, l/(4v^)), 
(0,4/^3) and (-1/4,1/(4^3)) (the green area in Fig. 1). A state is at most of W type 
(yellow in Fig. 1) iff its coordinates lie in the convex area enclosed by the lower edges 
of the triangle and the curve {xb{v), yB^v)) given by 

5(4 _ ^2) ' «-v-/ 4 4_^2 

with —1 < V < 1. All these boundary lines belong to the lower of the two 
entanglement classes which they separate (note that the fully separable state |(vrooo + 
TTiii) corresponding to (0, v^/4) is part of all three boundaries). 

An arbitrary state p can be projected onto the set of GHZ-symmetric states by the 
twirling operation 

S{p) = [ dU UpU^ (8) 



^b(^^) = — — , yB{v) = 7^ (7) 



where the integral is understood to include also averaging over the discrete symmetries, 
and dU is the Haar measure for the GHZ symmetry group. The coordinates of the 
symmetrized state are 

_ 1 , , 

y = ^ (^Pooo.ooo + Piii.iii - ^ ■ (9) 
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Since GHZ symmetry operations do not change the entanglement class, and mixing can 
result only in the same or a lower class, the entanglement class of S{p) is at most that 
of p. 

Now we determine the set of GHZ-symmetric states with positive partial transpose 
(PPT). Due to the permutation invariance of the matrix it does not matter which 
partition we use for the partial transpose. To identify the PPT states, we note that 
the only non-diagonal elements of p{x,y) are pooo,iii = Piii.ooo = x. A partial transpose 
obviously does not modify the diagonal elements, however it moves those non-diagonal 
elements along the anti-diagonal. Since the diagonal elements of a density matrix are 
always positive, the only eigenvalues which can be negative are those of the 2x2 
submatrix containing x, that is, 

f^"^' 1 \ ] (10) 

The two eigenvalues of this matrix are | — '^^V ^ they are both nonnegative 

exactly for the polygon of fully separable states. An alternative way to obtain this result 
would have been the criterion by Kay [17]. 



3. Optimal linear witnesses on subsets of states 

Linear entanglement witnesses are an important tool to detect entanglement [H El [6] , 
see also [7] and references therein. They are defined as Hermitean operators W whose 
expectation value tr (Wp) is positive for all unentangled states, but negative for some 
entangled states. If, for some state p, the expectation value of W is negative the witness 
is said to detect p. 

There exists a straightforward extension of this concept to class- specific 
entanglement [3]: The expectation value of a class-specific witness is positive for any 
state not belonging to the corresponding entanglement class or classes, and negative 
for some states in those entanglement classes. For example, a witness for GHZ-type 
entanglement is negative only for states containing GHZ-type entanglement. 

A useful property of linear entanglement witnesses is that they can be specified, 
up to an irrelevant positive factor, solely by the set of states for which they vanish, 
i.e., a hyperplane in the space of density matrices. This is derived directly from the 
corresponding property of general Hermitean operators by noting that the expectation 
value of the completely mixed state must be positive. 

To simplify terminology, in the following we will refer to states of the entanglement 
type(s) which we want to detect as "interesting", whereas we will call the states which 
are unentangled or have the "wrong" type of entanglement "uninteresting" , respectively. 
In this spirit, the uninteresting states will always be those whose entanglement type is 
"too low" according to the hierarchy given in the introduction, and therefore the set of 
uninteresting states will always be convex. 
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An important goal is to find optimal entanglement witnesses, that is, witnesses 
which detect as many states as possible. Lewenstein et al. [6] defined a witness W to be 
optimal if there is no other witness which detects all the states detected by W and also 
some other states not detected by W. Of course, this definition applies to class-specific 
witnesses as well. 

However, an optimal witness as defined above is not necessarily optimal for a given 
set. This is because a certain state may be detected by different optimal witnesses. For 
exaniple, we will see that there are many optimal witnesses detecting the standard GHZ 
statqj, but almost all of them are not optimal for detecting the pseudo-pure states in 
equation ([3]). That is, there exist such mixtures which are of GHZ type, but which are 
not detected by those globally optimal witnesses. 

For this reason, we specifically study witnesses that are optimal for a given (convex) 
set of states. This concept has been investigated by Hulpke et al. in Ref. [I6]. To define 
such a witness, the definition by Lewenstein et al. |6j can be modified to consider only 
the states in that particular set. In the important case of mixtures of a given highly 
entangled state, such as GHZ, with some sort of noise described by another state, the 
relevant set is just a straight line in the space of density matrices (typically, with the 
entangled state in the interesting class, and the noise in the uninteresting class). In that 
case, this definition reduces to the simple criterion: A witness is optimal on a given line 
if there is no interesting state on that line for which the expectation value of the witness 
is positive (cf. also Fig. 2). 

Note that, using the definition above, a witness which is optimal for a given set 
needs not be globally optimal: If the "mixture line" crosses a corner in the boundary 
of uninteresting states, for that point there will be generally more than one witness 
fulfilling the restricted optimality condition. If, in addition, that state is also on the 
boundary of the set of all states, it may happen that not all of those witnesses are 
globally optimal (see Fig. 2). However, there is always a witness which is both optimal 
for a given set and globally optimal. 

Importantly, for any convex subset M' of the set M of all states, every witness 
W which is both optimal with respect to the subset M' and vanishes on at least one 
full-rank state in M' , is also globally optimal (i.e., with respect to M). This can be seen 
as follows: If the witness W were not globally optimal, there would be another witness 
W" that detects all the interesting states in M detected by W, and then some more. 
Let p be such a state with tr(Wp) > > tr(>V'p), i.e., p is not located at the detection 
border of W. Since W is optimal for M', both witnesses must detect exactly the same 
states in M'. Now consider a full-rank state a G M' so that tr(Wcr) = (which implies 
tr(VV"cr) = 0), and the set of density matrices p(A) = (1 — A)(t + \p with real parameter 
A. Since a is of full rank, p(A) is still a density matrix for sufficiently small negative A. 
However for A < 0, the state is detected by W but not by W, in contradiction to the 
assumption. 

§ We will refer to these witnesses (as introduced in [6]) as globally optimal, as opposed to the witnesses 
which are optimal only for a certain subset of all three-qubit mixed states. 




Figure 2. Illustration for different witnesses in a convex set M of states which is 
a union of the disjoint sets of "interesting" and "uninteresting" states (indicated by 
green and blue regions, respectively). Moreover, the circle represents a convex subset 
of states M' C M, and the right-lower border of the figure another convex subset 
M" C M. The orange line Wi shows a generic witness in M (not optimal with respect 
to any set) while the red line W2 stands for a witness which is both globally optimal 
in M and also optimal with respect to subset Af (but not with respect to M"). The 
brown line Wa is an example for a witness which is optimal with respect to the border 
M" but not globally optimal (and also not optimal with respect to the circle M'). 

We are specifically interested in witnesses that are optimal for the subset of GHZ- 
symmetric states. These states form a plane triangle that contains the completely 
mixed state. Hence, for any witness W detecting a GHZ-symmetric state, the set of 
GHZ-symmetric states with tr(>Vp) = is a straight line. As we will see, all those 
witnesses vanish on at least one full-rank GHZ-symmetric state, and are therefore also 
globally optimal. 

4. GHZ-symmetric witnesses 

We now consider a special set of witnesses, the GHZ-symmetric witnesses, which are 
defined by the condition that for any GHZ symmetry operation U, we have UWW = W. 
It is easy to see that such witnesses must have the form 

W = a]l + 6|GHZ+)(GHZ+| +c|GHZ_)(GHZ_| (11) 
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The set of states for which the GHZ-symmetric witness has zero expectation value is 
then computed as 



The reverse is also true: For any straight line which does not cross the set of uninteresting 
states, there exists a corresponding GHZ-symmetric witness which is unique up to a 
positive factor. Proof: From equation f[T^ . it is immediately obvious that for each line, 
there exists (up to a real factor) exactly one GHZ-symmetric Hermitean operator whose 
expectation value vanishes on that line. The condition that for the fully separable 
state Pmixed at the origin the witness has positive expectation value fixes the sign 
via the inequality a + (6 + c)/8 > 0. Next we show that any such line which does 
not cross the set of uninteresting states corresponds to an entanglement witness. To 
see this, consider a state p which is detected by the GHZ-symmetric operator W, 
i.e., tr(>Vp) < 0. Then the symmetrized state S{p) is also detected by W, since 



tr(>V5(p)) = JdU tT{WUpW) = JdU tT{WWUp) = J dU tr(>Vp) = tr(Wp). Now 



S{p) is an interesting state because it lies on the side of the line which is opposite 
to the completely mixed state, and the line by assumption does not cross the set of 
uninteresting states. As the entanglement class of p is at least that of S{p), the state p 
must also be an interesting state, which completes our proof. 

From this one-to-one correspondence between the straight lines considered in the 
proof above and GHZ-symmetric witnesses one also sees that for any witness detecting 
GHZ-symmetric states, there is (up to a positive factor) a unique GHZ-symmetric 
witness which detects exactly the same GHZ-symmetric states. Therefore, in order 
to find witnesses that are optimal for GHZ-symmetric states, it is justified to consider 
only GHZ-symmetric witnesses. From Fig. 1 it is readily seen that a witness is optimal 
for GHZ-symmetric states iff it is tangential to the border of uninteresting states. This 
is sufficient knowledge for the explicit construction of witnesses which are optimal for 
GHZ-symmetric states. 

5. Construction of optimal GHZ-symmetric witnesses 

Since the set of GHZ-symmetric states is symmetric under the local unitary operation 
crf^ corresponding to p{x,y) i— p{—x,y), it suffices to consider witnesses which detect 
states with x > (since all the states p(x = 0, y) are fully separable, there exist no 
witnesses detecting states with both signs of x). The corresponding witness for a; < 
can be constructed by simply exchanging b and c in equation ffTTj) . 

We start with the detection of states which are not separable, that is, the interesting 
states are those containing any type of entanglement, while the uninteresting states are 
the fully separable ones. Here the border is the straight line through the points (1/8, 0) 
and (0, v^/4), which has the equation 




(12) 



— X — 




(13) 
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From this, together with equation (fT^ . the unique optimal witness (with respect to the 
GHZ-symmetric states) 

Wbisep\sep = 11-4 |GHZ+)(GHZ+| + 2 |GHZ_)(GHZ_| (14) 

is derived. 

Next we consider the detection of tripartite entanglement, i.e., the interesting 
states are GHZ and W states, while the uninteresting states are the biseparable and 
fully separable ones. We mention that this case was studied extensively for various 
multipartite systems in Refs. [HI [I9] . The border between those two is again a straight 
line, with the equation 

3 

- 2x - VSy + - = . (15) 
In analogy with the preceding case, we obtain the unique optimal witness 

>VH/\bisep = \^- |GHZ+)(GHZ+| (16) 

This witness was actually found a long time ago |3]. It is known to be globally optimal. 

Finally, we study witnesses detecting GHZ-type states. Here the boundary is given 
by equation ([7]), which is a convex curve. Thus, each tangent to this curve corresponds 
to a witness which is optimal for mixtures of the GHZ state and a state on the axis 
a; = (and, since those mixtures have full rank, also globally for all three-qubit states). 

Due to the finite curvature of that boundary, the parameters of these witnesses 
depend on the point {x-b^Vb) where the "mixing line" intersects the boundary. The 
tangent to curve ([7]) at parameter value f = f o is given by the equation 

y'B{vo)x - x'B{yo)y + (i/B(^^o)a;B(^o) - XB{yo)y'B{yo)) = . (17) 

From this, we obtain for the coefficients of the optimal witness 

A ( (yB{vo) + x'^ivo) - XBivo)yBivo)^ (18) 

(19) 



4^3, 

1 1 \ 

^y'sM - ^x'jiivo)] 



c = A (^-^y'^ivo) - ^^b(^o)) (20) 

where A has to be chosen such that a + (6 + c)/8 > 0. Inserting ([7]), together with an 
appropriate choice of A gives 

Wghz\h.(^o) = l^-^2_ + 4 ^GHZ, - ^ + 4 • (2^^ 

For the important case of the standard GHZ state mixed with white noise, equation 
the mixing line crosses the boundary at Vq = fws ~ 0.980701 and thus the optimal 
witness in that case reads 

y^Guz\w{vws) ~ 0.75 11 - 0.999876 7rGHZ+ - 0.433327 ttghz. (22) 
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It is instructive to compare this with the projection witness 

>Vp,oj = ^]l-|GHZ+)(GHZ+| . (23) 

This witness, which happens to be GHZ diagonal as well, is derived from the W- 

type state with maximum overlap to GHZ, which is \W+ ) = (|+ H — ) + |H h) + 

+ (here |±) = (|0) ± \l))/V2). The image of the symmetrized state 

Pmax = S{7ry[r^ ) is located at the point (3/8, 1/(2a/3)) in the triangle ofFig. 1 where the 

GWL/W boundary reaches the lower-right border. This state has vanishing expectation 
value for the projection witness (123|) . therefore the projection witness is optimal for 
the lower-right border of the triangle, that is, for mixtures of GHZ with p^. However, 
it is not optimal for the entire set of GHZ-symmetric states, as can easily be seen by 
calculating the optimal witness from equation fl2T|) . Since at the right border, vq = 1, 
the optimal witness for all GHZ-symmetric states is 

>VGHz\iy(l) = ^11 - |GHZ+)(GHZ+| - ^ |GHZ_)(GHZ_| . (24) 
This witness detects more GHZ-symmetric GHZ-type states than the projection witness 

6. Conclusions 

In this article, we have described how the two-parameter family of three-qubit GHZ- 
symmetric states — whose SLOCC classes of entanglement are known exactly — can be 
used to derive optimal witnesses for the SLOCC classes of general three-qubit states. 
To this end, we have studied in detail the concept of optimal witnesses with respect to 
a certain subset of states, introduced in |T6]. We have derived analytically the solutions 
for a bisperability witness, a witness of genuine three-qubit entanglement and a family of 
optimal witnesses for GHZ-type entanglement. Our findings confirm (witness for genuine 
entanglement) and improve (detection of GHZ-type vs. H^-type class) the results for the 
well-known three-qubit entanglement witnesses. In particular, our approach produces 
a simple one-parameter optimization scheme for witnessing GHZ-type entanglement in 
noisy GHZ states that can be applied for the assessment of three-qubit entanglement in 
experimentally measured density matrices. 
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